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ABSTRACT 

We report observations of 6 cm, 3.6 cm, 1.3 cm, and 7mm radio continuum, 
conducted with the Very Large Array towards IRAS 18566+0408, one of the 
few sources known to harbor H 2 CO 6 cm maser emission. Our observations 
reveal that the emission is dominated by an ionized jet at cm wavelengths. 
Spitzer/IRAC images from GLIMPSE support this interpretation, given the pres- 
ence of 4.5 /iin excess emission at approximately the same orientation as the cm 
continuum. The 7 mm emission is dominated by thermal dust from a flattened 
structure almost perpendicular to the ionized jet, thus, the 7 mm emission ap- 
pears to trace a torus associated with a young massive stellar object. The H2CO 
6 cm maser is coincident with the center of the torus-like structure. Our observa- 
tions rule out radiative pumping via radio continuum as the excitation mechanism 
for the H 2 CO 6 cm maser in IRAS 18566+0408. 
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1. Introduction 

How massive stars (M > 10 M Q ) form is among the most important unsolved problems in 
astrophysics. Large observational and theoretical efforts in the last ten years are supporting 
the hypothesis that massive stars form via accretion of material through disks; however 
other possible mechanisms such as spherical accretion and coalescence cannot be discarded, 
mainly in the case of stars more massive than 20 M (e.g., Beuther et al. 2007, Cesaroni 
et al. 2007, Beltran et al. 2006). Evidence for massive star formation via accretion disks 
is mainly indirect, e.g., through the detection of jets and bipolar molecular outflows (e.g., 
Arce et al. 2007; Garay et al. 2003, 2007); however more direct evidence, such as flattened 
rotating structures around young massive stellar objects, has been observed toward a number 
of objects (for example, see the review by Cesaroni et al. 2007). In some cases, CH3OH, 
H 2 0, SiO and OH masers have also been reported to trace circumstellar disks around young 
massive stellar objects (e.g., see reviews by Cesaroni et al. 2007 and Fish 2007). 

Although the formation mechanism of massive stars is still not well understood, the 
sites where massive stars are being (or have been recently) formed are readily identifiable 
by the effects that young massive stellar objects have on their environments. Among such 
effects are massive molecular outflows, the formation of H II regions, hot molecular cores, as 
well as the presence of maser emission from several molecular transitions (e.g., Churchwell 
2002). The 6 cm K-doublet line of formaldehyde (H 2 CO) is one of the molecular transitions 
that has been detected as maser emission toward massive star forming regions. H2CO 6 cm 
masers have several characteristics that distinguish them from other astrophysical masers: 
(1.) they have been found toward a small number of regions (e.g., Araya et al. 2007b, 2004; 
Mehringer et al. 1995; Pratap et al. 1994; Forster et al. 1985), (2.) they are weaker than 
other astrophysical masers such as H 2 0, OH, and CH 3 OH masers (e.g., Araya et al. 2007b, 
2007c; Hoffman et al. 2003, 2007), and (3.) H 2 CO masers have been only detected close 
to very young massive stellar objects (such as hot molecular cores and deeply embedded 
infrared sources, e.g., Araya et al. 2006a) and appear not to be associated with high velocity 
outflows or more evolved phases of massive star formation such as compact H II regions (see 
review article by Araya et al. 2007a). 

As pointed out by Mehringer et al. (1995), the paucity of H 2 CO 6 cm masers suggests 
that the pumping mechanism of the masers may require quite specific physical conditions for 
the inversion, conditions that may be short lived in massive star formation environments (see 
Araya et al. 2007a for a discussion of other reasons why H 2 CO masers may be rare). If this 
were the case, then H 2 CO 6 cm masers could become a useful probe for detection of specific 
physical conditions; however at present there exists no successful model to explain H 2 CO 
6 cm masers. Although infrared and collisional excitation mechanisms have been discussed in 
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the literature as agents capable of inverting the H2CO 6 cm levels (see Litvak 1970, Thaddeus 
1972), only Boland & de Jong (1981) developed a model specifically intended to explain one 
of the known H 2 CO 6 cm maser sources, i.e., the maser in NGC 7538. Their model is based on 
radiative pumping of the H 2 CO 6 cm levels by radio continuum from a background compact 
(emission measure > 10 8 pccm~ 6 ) H II region. 

Araya et al. (2005) explored whether the Boland Sz de Jong (1981) model is capable 
of explaining the H 2 CO 6 cm maser in IRAS 18566+0408, but could not draw a strong con- 
clusion due to insufficient observational constraints on the radio continuum properties of 
IRAS 18566+0408. Thus, to more critically test whether a pumping mechanism by radio 
continuum can explain the H 2 CO maser in IRAS 18566+0408, and to explore the relation 
between the H 2 CO maser and the massive star formation process in the region, we conducted 
a multi-wavelength study of the radio continuum in IRAS 18566+0408. In §1.1 we summa- 
rize the properties of IRAS 18566+0408, our observations and data reduction procedure are 
reported in §2; the results and discussion of the radio continuum observations are presented 
in §3. In §4 we report our conclusions. 



1.1. IRAS 18566+0408: A Massive Disk Candidate with an H 2 CO Maser 

Located at a distance of 6.7 kpc, IRAS 18566+0408 (G37.55+0.20; Mol 83, Molinari et 
al. 1996) is a massive star forming region with a bolometric luminosity of ~ 6 x 10 4 Lq, 
equivalent to the luminosity of an 08 ZAMS star (Sridharan et al. 2002; Araya et al. 2005)|j. 
Current star formation in the region is evident from the presence of multiple molecular 
outflows, as well as 22 GHz H 2 and 6.7 GHz CH 3 OH maser emission (Beuther et al. 2002a, 
2002b). Zhang et al. (2007) found that roughly all the far IR luminosity comes from a single 
compact (< 5") source, indicating the presence of an embedded massive protostar. 

IRAS 18566+0408 is an ideal source for the study of massive star formation because 
only a single weak radio continuum source is found in interferometric observations (Carral et 
al. 1999, Araya et al. 2005), and thus the study of the region is less susceptible to confusion. 
The absence of strong radio-continuum indicates that IRAS 18566+0408 is in a phase prior 
to the development of a bright ultra-compact H II region. IRAS 18566+0408 also harbors 
one of the few known H 2 CO 6 cm masers in the Galaxy (Araya et al. 2005). 

Zhang (2005) listed IRAS 18566+0408 as a massive circumstellar disk candidate. Ac- 



1 Recently, Zhang et al. (2007) estimated a total luminosity of ~ 8 x 10 4 L Q mainly based on mid to far 
IR data at ~ 1' resolution. 
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cording to Zhang (2005), the massive dislo has a physical size of 8000 AU and mass of 60 M Q . 
Zhang et al. (2007) have recently published a multi-frequency study of IRAS 18566+0408. 
In §3.3 we discuss their results in the context of the observations reported in this paper. 

2. Observations and Data Reduction 

We report VLA§| C (6 cm), X (3.6 cm), K (1.3 cm), and Q (7mm) band radio continuum 
observations of IRAS 18566+0408. The C and X band observations were carried out on 
2005 May 18, with the VLA in the B configuration. The K and Q band observations were 
conducted on 2005 August 5 and 2003 April 25 with the VLA in the C and D configurations, 
respectively. The default VLA continuum mode (4IF, 50 MHz per IF) was used for the 
observations. 

The total integration time on-source was approximately 1.6, 1.0, 1.9, and 0.5 hours 
for the C, X, K, and Q band observations, respectively. The observations were conducted 
in phase referencing mode, with switching cycles (source/calibrator) of 720/120, 600/150, 
120/40, and 50/50 seconds for the C, X, K, and Q band observations, respectively. We 
used the fast-switching VLA observing mode for the K and Q band observations, given 
that short calibration cycles are required to track rapid phase variations due to tropospheric 
instabilities. 

The data reduction was conducted following the VLA high-frequency recipe for calibra- 
tion and standard imaging procedures using the NRAO package AIPS. In Table 1 we give 
the details of the observations. 



3. Results & Discussion 

We detected radio continuum emission from IRAS 18566+0408 at all observed wave- 
lengths. In Figure 1 we show the images obtained in the different bands. In the figure 
we also mark the position of the H 2 CO 6 cm maser (RA = 18 h 59 m 09 f 5975±0 f !002, Decl. = 



2 In this paper we use the term "massive disk" as synonymous of "circumstellar torus" (also referred as 
pseudo-disks in the literature, Cesaroni et al. 2007). The term "massive disk" should not be confused with 
"accretion disk" . 

3 The Very Large Array (VLA) is operated by the National Radio Astronomy Observatory (NRAO), a 
facility of the National Science Foundation operated under cooperative agreement by Associated Universities, 
Inc. 
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+04 o 12 / 15"57±0 , . / 03, J2000, Araya et al. 2005). In Table 2 we list the details of the radio 
images and the parameters of the continuum emission in IRAS 18566+0408. 

3.1. An H 2 CO 6 cm Maser Pinpointing a Possible Circumstellar Torus 

In Figure 2 we show the radio spectral energy distribution (SED) of IRAS 18566+0408 
based on our observations plus the 2 cm detection by Araya et al. (2005). We found that the 
SED is satisfactorily fitted by the combination of optically thin thermal dust emission plus 
optically thin Bremsstrahlung radiation from an ionized jet. The dust emission in Figure 2 
was fitted assuming a dust temperature of 50 K, i.e., approximately the temperature of 
the warm dust component reported by Sridharan et al. (2002), and a dust opacity k v 
proportional to A~ L8 (with K43.3GHZ — 0.044 cm 2 g _1 ) which was obtained from a fit of the 
opacities for dust particles with thin ice mantles and n# 2 = 10 6 cm" 3 presented by Ossenkopf 
& Henning (1994). Our assumption of k v proportional to X' 13 where (3 = 1.8, is consistent 
with the recent empirical determination of (3 by Hill et al. (2006) for dust at 50 K. 

Although the cm emission can be also well fit by free-free radiation from a simple 
spherical HII region (see §3.2), we decided to fit the cm emission using an ionized jet model 
given the elongated morphology of the 6 cm and 3.6 cm continuum emission (Figure 1); 
specifically, we used the "standard" collimated model of Reynolds (1986). The fit was 
obtained assuming an electron density of 2xl0 4 cm -3 at 1340 AU (i.e., 0'.'2 at 6.7kpc) from 
the protostar, an electron temperature of 10 4 K, fully ionized gas, a width at the base of 
the jet of 1340 AU, and a half-length of the ionized jet of 6700 AU. In addition, we assumed 
that the jet is in the plane of the sky. The estimated ionized mass-loss rate is greater than 
10~ 6 M Q yr _1 (assuming a jet velocity of >100 km s _1 ), which, for example, is comparable to 
the mass-loss rate of the ionized jet in GGD 27 (Marti et al. 1995). The overall cm- mm SED 
of IRAS 18566+0408 is very similar in shape, but much more luminous, when compared to 
those found in some low mass young stars (e.g., L1551 IRS5, Rodriguez et al. 1998). Since 
we only detected optically-thin radio emission from the jet, we cannot discriminate between 
different Reynold's models, i.e., by adjusting the parameters all models can produce good 
fits. Thus, we are unable to precisely determine physical parameters of the ionized gas such 
as the electron temperature and density. 

The interpretation that the elongation detected at 6 cm and 3.6 cm is caused by an 
ionized jet is strengthened by Spitzer/IRAC data of IRAS 18566+0408 from GLIMPSE. In 
Figure 3 we show a three-color image of IRAS 18566+0408 in the 3.6 /xm (blue), 4.5 /zm 
(green), and 8.0 /im (red) bands. The Spitzer data show an excess at 4.5 ^m toward the 
north-west of the IRAS 18566+0408 position. Excess in the 4.5 /iin band is a known tracer 
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of outflows (shocked gas) in massive star forming regions due to the contribution of H 2 
emission lines (e.g., Noriega-Crespo et al. 2004, Smith et al. 2006). Thus, the similar 
orientation of the 4.5 /im excess and the cm radio continuum supports the interpretation 
that the elongation observed at 6 cm and 3.6 cm is due to an ionized jet. 

In Figure 4 we show a three-color image of the 6 cm (blue), 3.6cm (green), and 7mm 
(red) emission as shown in Figure 1. The ionized jet is elongated in a ~E-W orientation, 
and although the signal-to-noise ratio of the Q-band detection is low (~ 7er), the orientation 
of the 7 mm source (after beam deconvolution) is in the NE-SW direction. Since the 7 mm 
emission is almost perpendicular to the ionized jet and is dominated by thermal dust, it 
appears that the 7mm source traces a circumstellar torus in IRAS 18566+0408. In Figure 4 
we also show the location of the H 2 CO 6 cm maser (Araya et al. 2005). The maser is 
coincident with the center of the torus-like structure. 

Using the 7 mm dust flux density from the fit (Figure 2), assuming a 7 mm dust opacity 
of «43.3GHz = 0.044 cm 2 g _1 , a gas-to-dust ratio of 100, and a dust temperature of 50 K, the 
mass of the 7 mm source is ~5OOM (if the gas temperature is 100 K, the mass would be a 
factor of 2 smaller). Thus, the 7mm source is indeed tracing a massive flattened structure. 
This torus candidate appears to be more massive and greater than the parameters reported 
by Zhang (2005; see however §3.3). Our mass determination should be considered a rough 
estimate because of the uncertainties in the assumed gas-to-dust ratio, dust opacity, and 
dust temperature values. Nevertheless, we find that the mass of the 7 mm source must 
be larger than 10 M , since even after assuming an extremely large 7mm dust opacity 
(/«43.3GHz = 0.78cm 2 g _1 , Rodriguez et al. 2007), the mass is still > 10 M Q . Assuming that 
a single 08 ZAMS star is forming in IRAS 18566+0408 and a torus mass > 10 M , then 
the M star /M torus ratio for the system would be < 2, suggesting that the torus-like structure 
is dynamically unstable (e.g., see Rodriguez et al. 2007, and references therein). Further 
high angular resolution observations of the region are needed to investigate the stability and 
kinematics of the 7 mm source. 



3.2. Applicability of the Boland & de Jong (1981) Model 

To test whether population inversion by radio continuum radiation can explain the 
H 2 CO maser in IRAS 18566+0408, we adapted the Boland & de Jong (1981) formalism to 
the specific characteristics of IRAS 18566+0408. In particular, we included the contribution 
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of dust emission to the photon occupation numbecl and we used the following expression for 
the optical depth of the radio continuum emission (see Reynolds 1986): 

r v (y) = 2a k w n 2 x 2 T- 1 - 35 (y/y o y^~ 2 -\ (1) 

where a k = 0.212, q T = —2 for the "standard" collimated model of Reynolds (1986), w a = 
670 AU, n Q = 2 x 10 4 cm~ 3 , x Q — 1, T — 10 4 K (see §3.1). We assumed that the ionized jet 
is in the plane of the sky (i.e., sin(i) = 1 in the Reynolds 1986 notation, see also Figure 1 in 
Reynolds 1986), and hence y is simply the distance along the ionized jet, and y Q is the distance 
at which the flow is injected (we assumed y Q = 1340 AU based on the IRAS 18566+0408 cm 
emission images, Figure 1). 

Given the observational constraints imposed by our data, we find that the model does 
not predict population inversion of the 6 cm K-doublet level. The specific ionized gas model 
used to fit the cm radio continuum is not relevant to the outcome of the maser model as 
long as the free-free emission is optically thin at A < 6 cm as observed. For example, the 
cm SED of IRAS 18566+0408 can be fitted by optically thin emission from a cylindrical H II 
region characterized by n e = 1.7 x 10 3 cm -3 , T e = 10 4 K and EM= 1.4 x 10 5 pccm~ 6 . No 
inversion is predicted for such a background H II region (see discussion of Figure 1 in Boland 
& de Jong 1981, and Figure 4 in Araya et al. 2005). We conclude that excitation by radio 
continuum is not a feasible mechanism to explain the maser in IRAS 18566+0408. 

3.3. Comparison with the Observations by Zhang et al. (2007) 

Recently, Zhang et al. (2007) reported a multi-wavelength study of IRAS 18566+0408. 
They conducted 22, 43, and 87 GHz radio continuum observations, as well as NH 3 (1,1), 
(2,2), and (3,3), SiO (2-1), and HCN (1-0) observations with the VLA and OVRO. In this 
section we discuss the consistency of the results found in our work with respect to the results 
found by Zhang et al. (2007). 



4 In the case of NGC 7538, Boland & de Jong (1981) neglected dust emission because they considered that 
free- free radio continuum radiation was dominant at all wavelengths of interest. 
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3.3.1. A Jet-like Outflow in IRAS 18566+ 04 08 

The SiO (2-1) observations by Zhang et al. (2007) reveal a well collimated molecu- 
lar outflow in IRAS 18566+0408. The outflow appears to be emanating from the torus-like 
structure reported in this work. The SiO outflow is in the SE-NW direction which is consis- 
tent with the orientation of the ionized jet and the large structure outflow traced by Spitzer 
(see Figures 3 and 4). However, the SiO outflow is not strictly collinear to the ionized jet 
or the 4.5 \xm excess emission: within ~ 1" from the center of the torus-like structure (Fig- 
ure 4), the ionized jet appears to have a position angle of ~ 100°; at ~ 10" from the torus, 
the position angle of the SiO outflow is 135°; and at ~ V scales, the 4.5 /xm outflow has a 
position angle of ~ 110°. Given that the velocity gradient of the SiO outflow is opposite of 
that of the CO outflow, Zhang et al. (2007) mentioned that the outflow may be very close 
to the plane of the sky and that the velocity inversion may be due to precession. The change 
in position angle between the radio jet, the SiO and the 4.5 \im outflow may also be due to 
precession of the jet. 

The NH3 and HCN data reported by Zhang et al. (2007) show an elongation parallel 
to the SiO outflow, indicating that these molecules may also be tracing the outflow. The 
NH 3 (1,1) and (2,2) velocity integrated emission is distributed almost symmetrically with 
respect to the center of the torus candidate, whereas the SiO, NH 3 (3,3) and HCN integrated 
emission is elongated mostly toward the NW of the 7 mm source. The asymmetry of the 
SiO, NH 3 (3,3) and HCN emission with respect to the position of the torus-like structure 
is similar to the one shown by the Spitzer 4.5 [im excess, in which 4.5 fim excess emission 
is only detected toward the NW of the radio continuum source. Density gradients in the 
medium may be responsible for such an asymmetry. 

3.3.2. Continuum Observations 

Zhang et al. (2007) conducted continuum observations at 22, 43, and 87 GHz. At 87 GHz 
{0 S yn ~ 5") they detected a continuum source coincident with the weak radio continuum 
source reported by Carral et al. (1999) and Araya et al. (2005, see also Figure 1). The 
peak intensity of the emission is 18mJybeam _1 and the total flux density of the source is 
31 mJy The 87 GHz flux is clearly dominated by dust emission, which is consistent with 
our radio SED interpretation. Zhang et al. (2007) detected a second mm source (MM-2) 
north-west from the position of the protostellar candidate (see their figure 1). We did not 
detect a counterpart of MM-2 in our observations, which is consistent with the non-detection 
at 43 GHz by Zhang et al. (2007). 
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The 43 GHz observations reported by Zhang et al. (2007) were conducted with the VLA 
in the DnC array, resulting on a synthesized beam of 2'.'7 x l'/3, which is more elongated than 
our 9 syn at 43 GHz (Table 2). The observations of Zhang et al. (2007) have a smaller rms than 
that of our observations, i.e., 0.1 versus 0.17mJy beam -1 . As in our observations, they also 
detected 43 GHz emission coincident with the radio continuum source. The source is slightly 
resolved with a peak intensity of 1 mJy beam -1 and flux density of 1.7 mJy The deconvolved 
size of the emission is 2'/0 x 1'.'2, P.A. = 24°, which corresponds to a major axis of 1.3 x 10 4 AU. 
Thus, the 7 mm source detected by Zhang et al. (2007) appears to be more compact and 
weaker than the 7 mm emission reported in this work (Table 2). To investigate the origin of 
this discrepancy we reduced the 7 mm observations by Zhang et al. (2007) from the VLA 
archive, and found that the discrepancy is mainly due to different weighting schemes used in 
the imaging, i.e., our 7mm image was created with natural weighting to improve signal-to- 
noise, whereas Zhang et al. (2007) apparently used a more uniform weighting scheme given 
that they have better sensitivity data. For example, if our data are imaged using a Briggs' 
robust parameter of 0, then the peak intensity of the 7 mm source is 1.0 mJy beam -1 (i.e., the 
same value reported by Zhang et al. 2007) and the integrated flux density is 2.6±0.6mJy. 
Our reduction of the Zhang et al. data using a robust weighting resulted in an integrated 
flux density of 1.74±0.35 mJy, i.e., we reproduced the value reported by Zhang et al. (2007) 
and the flux density is consistent within the errors with the corresponding flux density 
derived from our data. Regarding the size and position angle, given the low signal-to-noise 
of both data sets, and that the 7 mm source is barely resolved, then it is difficult to accurately 
measure the deconvolved size and position angle of the source. However, for example, if a 
robust weighting is used in the imaging of the Zhang et al. (2007) data, we found that 
the maximum size of the 7 mm source is 3'/7 x l'/7 with a position angle of 44°, which is 
consistent with our results (Table 2). Higher sensitivity and angular resolution observations 
are needed to precisely measure the size and position angle of the 7 mm source. Nevertheless, 
our observations and those of Zhang et al. (2007) are consistent with an elongated source 
almost perpendicular to the ionized jet that we report in this work. 

Zhang et al. (2007) also conducted VLA continuum observations at 22 GHz and found 
no radio continuum emission at a la level of 0.14mJybeam _1 . They also pointed out that 
Miralles et al. (1994) did not detect 2 cm emission at a la level of 0.16 mJy whereas Araya et 
al. (2005) report a 0.7 mJy source at 2 cm. Zhang et al. (2007) mentioned that variability of 
the radio continuum emission in the source may be responsible for this apparent discrepancy. 
We note that Miralles et al. (1994) did not consider a la value as an appropriate detection 
limit for the radio continuum, but they report upper limits of < 0.5 and < 0.8 mJy for 
6 cm and 2 cm emission which are consistent with our 6 cm detection (Table 2) and the 2 cm 
detection by Araya et al. (2005). Regarding the non-detection of 22 GHz emission reported 
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by Zhang et al. (2007), the 5a detection limit of their data is 0.7mJybeam _1 , and the 
peak intensity of the 22 GHz source reported in this work (see Figure 1) is 0.65 mJy beam -1 , 
i.e., at a 5a detection limit, the two data sets are consistent. The non-linear nature of 
interferometric observations, particularly at K band were tropospheric instabilities may lead 
to phase decorrelation, can be responsible for the non-detection of the source at a 5a level. 
Thus, given the available radio continuum data, we find no strong evidence for variability of 
the source. Nevertheless, the occurrence of flares in the H 2 CO maser (Araya et al. 2007c, 
see also Araya et al. 2007a) may be caused by variability of the radio continuum, thus future 
observations are required to investigate whether, and up to what level, the radio continuum 
emission is variable in this source. 



3.3.3. Mass Determination 

Based on mm and submm wavelength data, Zhang et al. (2007) report a spectral index 
a of 3.9 (Sj, oc u a ) and thus an emissivity index (3 = 1.9, which is consistent with the 
value of P — 1.8 derived from theoretical as well as empirical considerations (see §3.1). 
However, in an effort to better constrain the value of (3 by matching beam sizes, Zhang et 
al. (2007) derived a value of (3 = 1.3 which they believe to be more appropriate for the 
dust mass determination. Based on NH 3 observations, they obtained a kinetic temperature 
of 80 K for the molecular core associated with the 7 mm, which was assumed to be the dust 
temperature. In addition they assumed an opacity /t(250/im) = 12cm 2 g~ 1 . Using these 
parameters, Zhang et al. (2007) estimated a mass of 70 M© within 30000 AU from the 7 mm 
source. Using the same parameters (T rf = 80 K, /t(250/im) = 12cm 2 g -1 , j3 = 1.3) we derive 
a mass of ~ 80 M based on our observations. Given the flux density uncertainties (Table 
2) our derived mass is consistent with the 70 M mass reported by Zhang et al. (2007). 

3.3.4- A Molecular Core coincident with the Torus Candidate 

Zhang et al. (2007) detected compact molecular emission from the position of the torus 
candidate in their NH 3 and HCN observations. As mentioned above, the NH 3 emission at the 
position of the 7 mm source traces warm (~ 80 K) material, and in addition, the line width 
of the NH 3 emission is substantially broader (FWHM = 8.7km s" 1 ) in comparison with the 
more extended NH 3 emission (FWHM < 2km s _1 ). Zhang et al. (2007) find it unlikely that 
the broad linewidth is due to the outflow, instead the large linewidth may be due to relative 
motion of multiple cores, or due to infall/rotation. They found that if rotation is assumed, 
then the dynamical mass (assuming gravitationally bound motion) would be similar to the 
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70 M mass reported by Zhang et al. (2007, see §3.3.3). However, their data do not show 
direct kinematic evidence for rotation of the core. 

The peak velocity of the NH 3 emission is 85.2 km s" 1 , which closely agrees with the 
systemic velocity of the source as traced by CS emission and H 2 CO absorption (Bronfman et 
al. 1996; Araya et al. 2004). The LSR velocity of the H 2 CO maser in the region is ~ 80 km 
s -1 , thus the maser is not tracing gas at the systemic velocity of the molecular core but is 
tracing blueshifted material, perhaps associated with the base of the molecular outflow or a 
flaring disk with radial velocity gradients. 

4. Summary 

We report 6 cm, 3.6 cm, 1.3 cm, and 7 mm observations of IRAS 18566+0408 conducted 
with the VLA. We detected radio continuum emission from the source at all wavelengths. 
The cm emission appears to trace an ionized jet, whereas the 7 mm band is dominated by 
thermal dust emission. The interpretation that the cm radiation arises from an ionized jet is 
strengthened by Spitzer/IRAC observations of IRAS 18566+0408, showing excess at 4.5 //m 
at the same orientation as the suggested ionized jet. 

The 7 mm structure is elongated almost perpendicular to the ionized jet. Taking into 
account that the 7 mm emission is dominated by thermal dust, then the 7 mm source appears 
to trace a circumstellar torus. The mass of the torus-like structure is greater than 10 M Q (and 
as large as several hundreds solar masses depending on the assumed dust opacity and tem- 
perature). Given that definitive identification of a disk/torus requires kinematical evidence 
of rotation, future sub-arcsecond angular resolution molecular line observations are required 
to confirm the torus interpretation presented in this paper. 

Our results are consistent with recent molecular and continuum observations reported by 
Zhang et al. (2007). They found a SiO outflow that is approximately parallel to the ionized 
jet reported in this work, and a clump of NH 3 and HCN molecular gas at the position of 
the torus-like structure. Zhang et al. (2007) also conducted 7 mm VLA observations of the 
region and found an elongated structure whose orientation is almost perpendicular to the 
ionized jet reported in this work, which confirms our 7 mm results. 

The H 2 CO 6 cm maser in IRAS 18566+0408 (Araya et al. 2005) is coincident with the 
center of a torus-like structure which harbors a massive protostellar candidate. This result 
strengthens the association of H2CO 6 cm masers with very young massive stellar objects. 
The observed radio spectral energy distribution rules out pumping of the H 2 CO 6 cm maser 
in IRAS 18566+0408 by radio continuum, which Boland & de Jong (1981) proposed as the 
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excitation mechanism for the H 2 CO maser in NGC 7538. 
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Table 1. VLA Observations 



Parameter 


C-Band (6 cm) 


X-Band (3.6 cm) 


K-Band (1.3 cm) 


Q-Band (7 mm) 


Date 


2005 May 18 


2005 May 18 


2005 August 5 


2003 April 25 


VLA Configuration 


B- Array 


B- Array 


C- Array 


D- Array 


vo (GHz) 


4.86 


8.46 


22.46 


43.34 


RA a 


18 h 59 m 09?9 


18 h 59 m 09?9 


18 h 59 m 10?0 


18 h 59 m 09?9 


Dec a 


04° 12' 10" 


04° 12' 10" 


04°12'14" 


04°12'14" 


Flux Density Calib. 


3C286 


3C286 


3C286 


3C286 


Assumed (Jy) 


7.49 


5.21 


2.52 


1.45 


Phase Calib. 


J1824+107 


J1824+107 


J1851+005 


J1851+005 


Measured S„ (Jy) 


0.77 


0.72 


0.97 


0.66 



a Phasc tracking center (J2000). 
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Table 2. Parameters of the Radio Continuum Images 



Parameter 


C-Band (6 cm) 


X-Band (3.6 cm) 


K-Band (1.3 cm) 


Q-Band (7 mm) 


Syn. Beam 


1'.'3 x l'/2 


0'.'93 x OT/79 


l'/4 x l'.'O 


2"0 x 1'.'6 


Syn. Beam P.A. 


-27° 


24° 


5° 


22° 


rms (/iJy b -1 ) 


24 


21 


93 


168 


(mjy) a 


0.52±0.09 


0.55±0.10 


0.80±0.18 


3.1±0.6 


Peak RA a 


18 h 59 m 09?97±0?01 


18 h 59 m 09?99±0?01 


18 h 59 m 09?98±0.01 


18 h 59 m 09?99±0?01 


Peak Dccl. a 


04°12'15'.'49±tf'08 


04°12'15'.'67±tf'08 


04°12'15'.'74±0.06 


04°12'15'.'7±0'/2 


Deconvolved Size" 


2'.' 9 x (X'73 


2'.'3 x 0{'85 


< 1"5 x 1'.'5 


3{'3 x 1'.'6 


Deconvolved Size P.A.° 


85° 


103° 




45° 


Physical Size (10 3 AU) 


19 x 5 


16 x 6 


< 10 x 10 


22 x 11 



a Parameters obtained from a 2D Gaussian fit of the brightness distribution using the task JMFIT in AIPS. 



-17- 



04 1 2 20 
19 
18 

I " 

; 16 
t 

: 15 



o 



10.2 



6 cm Cont. 




10.1 10.0 09.9 09.8 09.7 
RIGHT ASCENSION (J2000) 
Peak Intensity = 0.19mjy/b 
Levs = 0.024 X (-3, -2.5, 2.5, 3, 4, 5, 6, 7, 7.8) mjy/b 



04 12 20 

ig - 

18 - 



3.6 cm Cont.j^ 



I H I I I I L_ 

59 10.3 10.2 10.1 10.0 09.9 09.8 
RIGHT ASCENSION (J2000) 

Peak Intensity = 0.17mjy/b 

Levs = 0.021 X (-3, -2.5, 2.5, 3, 4, 5, 6, 7, 7.9) mjy/b 



1 .3 cm Cont. 

o 




<3 



39 10.3 1 0.2 1 0.1 10.0 09.9 09.S 
RIGHT ASCENSION (J2000) 
Peak Intensity = 0.65 mjy/b 
Leys = 0.091 x (-3, 3, 4, 5, 6, 6.9) mjy/b 



04 12 20 

19 - 
18 - 



10.2 



7 mm Cont. 




O 



o 



10.1 10.0 09.9 09.8 
RIGHT ASCENSION (J2000) 
Peak Intensity = 1.25 mjy/b 
Levs = 0.17 X (-3, -2.5, 2.5, 3, 4, 5, 6, 7, 7.4) mjy/b 



Fig. 1. — Radio continuum emission detected toward IRAS 18566+0408 at the four different 
wavelengths observed with the VLA. In the lower left corner of each image we show the 
respective synthesized beam (the size and position angle of the synthesized beams are given 
in Table 2). The location of the H 2 CO 6 cm maser imaged with the VLA-A by Araya et al. 
(2005) is indicated by "©". The size and position angle of the "©" symbol represents the 
size and position angle of the synthesized beam of the Araya et al. (2005) observations. 
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Fig. 2. — Radio SED of IRAS 18566+0408 based on the observations reported in this paper 
(Table 2), in addition to the 2 cm data point by Araya et al. (2005). Dashed line is fit 
to the ionized gas emission using the "standard" collimated ionized jet model developed by 
Reynolds (1986). At millimeter wavelengths the SED is dominated by optically thin thermal 
dust emission. The total fit to the SED, i.e., the addition of the ionized gas and thermal 
dust flux density, is shown with a continuous line. 
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Fig. 3. — Three-color image of IRAS 18566+0408 from the Spitzer/IRAC 3.6 /im 
(blue), 4.5/zm (green), and 8.0/im (red) bands. The IRAC detectors were saturated 
at the IRAS 18566+0408 peak. Note the 4.5 fim (green) excess toward the NW of 
IRAS 18566+0408. 4.5 fim excess is a known tracer of shocked gas in massive star form- 
ing regions (e.g., Smith et al. 2006). In the inset we show the 3.6cm radio image obtained 
here (see Figure 1). The elongation of the 3.6 cm continuum source is approximately parallel 
to the 4.5 /iin NW excess, which agrees with the suggestion that the ionized gas traces a jet 
in IRAS 18566+0408. 
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Fig. 4. — Three-color image of the 6 cm (blue), 3.6 cm (green), and 7mm (red) radio contin- 
uum emission from IRAS 18566+0408. The colors were saturated to increase contrast. At 
cm wavelengths (where the emission is dominated by ionized gas) the brightness distribution 
is oriented East- West (highlighted by thick white lines), whereas at 7mm the orientation is 
NE-SW. Given that the 7 mm emission is dominated by thermal dust and the ionized jet 
is almost perpendicular to the 7 mm source, then the 7 mm flattened structure appears to 
trace a circumstellar torus. The H 2 CO maser position as detected by Araya et al. (2005) is 
shown by "©". The maser is coincident with the center of the 7 mm source. 



